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Relativistic Wavepackets in Classically Chaotic Quantum Cosmological Billiards 
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Close to a spacelike singularity, pure gravity and supergravity in four to eleven spacetime di- 
mensions admit a cosmological billiard description based on hyperbolic Kac-Moody groups. We 
investigate the quantum cosmological billiards of relativistic wavepackets towards the singular- 
ity, employing flat and hyperbolic space descriptions for the quantum billiards. We find that the 
strongly chaotic classical billiard motion of four-dimensional pure gravity corresponds to a spreading 
wavepacket subject to successive redshifts and tending to zero as the singularity is approached. We 
discuss the possible implications of these results in the context of singularity resolution and compare 
them with those of known semiclassical approaches. As an aside, we obtain exact solutions for the 
one-dimensional relativistic quantum billiards with moving walls. 
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I. INTRODUCTION 

Spacetime singularities are known to generically ap- 
pear in the classical theory of general relativity [T] , and 
they are expected to be resolved through quantum ef- 
fects. Understanding the classical structure near these 
singularities should shed light on the mechanisms that 
lead to the quantum resolution. Via the Belinskii- 
Khalatnikov-Lifshitz (BKL) analysis [HIS] (see also 0J), 
the equations of general relativity can be traced back to- 
wards the initial singularity, leading in the limit to an 
ultralocal description of spacetime in terms of a system 
of steep walls with the logarithms of the spatial scale 
factors being effectively described by a massless particle 
bouncing and specularly reflecting off these walls. In be- 
tween these reflections, the universe is undergoing free 
Kasner evolutions [SJ [B] . Generally, a dynamical system 
of a pointlike particle moving along geodesies on a Rie- 
mannian manifold and bouncing specularly off its piece- 
wise smooth boundary is called a billiard system, and 
thus the relativistic billiards arising from the BKL ana- 
lysis have been termed cosmological billiards. Depending 
on the possible matter content and the dimension of the 
respective gravity theory, this motion may or may not 
be chaotic [Zl [8]. If it is chaotic, as is the case e.g. for 
pure gravity in four and supergravity in eleven spacetime 
dimensions, then the motion in scale factor space under- 
goes an inflnite series of epochs of anisotropic expansion 
and contraction as the singularity is approached. 

While the dynamics of the gravitational field generally 
shows a hidden symmetry in terms of Lorentzian Kac- 
Moody algebras [5], an underlying algebraic structure of 
pure gravity and supergravity theories with a number of 
spacetime dimensions from four to eleven was also re- 
vealed through the shapes of the billiard tables which 
arise in the BKL limit. Namely, these billiard tables co- 
incide with the fundamental Weyl chambers of infinite- 
dimensional Kac-Moody algebras [71 [5] , the latter being 
of the hyperbolic type if the corresponding billiard mo- 
tion is chaotic, as is the case for all the maximally su- 



persymmetric theories that are candidates for a unified 
description of the fundamental forces. Remarkably, the 
highest possible rank for a hyperbolic Kac-Moody alge- 
bra implies a maximum of eleven spacetime dimensions 
for the corresponding supergravity theory. 

For theories in a certain number of spacetime dimen- 
sions, the even Weyl groups of the respective algebra have 
in turn been identified with generalized modular groups 
based on arithmetic integer rings of different kinds of di- 
vision algebras [TOl[Tl]. For example, the billiard domain 
corresponding to Z) = 4 pure gravity is given by the fun- 
damental Weyl chamber of the infinite-dimensional hy- 
perbolic Kac-Moody algebra AE^ [T2], or equivalently 
by half the fundamental domain of the standard mod- 
ular group on the hyperbolic upper half-plane. In the 
case oi D = 11 supergravity on the other hand, the rele- 
vant Kac-Moody algebra is Eiq [9" , and the billiard table 
is given on the nine-dimensional generalized upper half- 
plane by half the fundamental domain of the generalized 
modular group with respect to integer octonions, see [T3] 
for details on the explicit construction of such billiard 
domains. 

For the resolution of the initial singularity of the classi- 
cal gravity theory, quantum effects are expected to be the 
essential additional ingredient. First insights in this di- 
rection arise from the quantization of the relativistic cos- 
mological billiards system. Following the original sugges- 
tion of 55 , the potential walls are translated into Dirich- 
let boundary conditions on the wavefunction, leading to 
a quantum cosmological billiard description evolving ac- 
cording to the free Wheeler-DeWitt equation and sub- 
jected to the boundary conditions. The generally con- 
stantly moving domain walls can be turned into static 
ones through a transformation to hyperbolic space, lead- 
ing to an eigenvalue problem of the Laplacian on the 
arithmetic billiard domain (see e.g. [15j for the math- 
ematical background in the case of complex numbers). 
While an exact solution of this classically chaotic quan- 
tum billiard system is unknown and out of reach, we 
pointed out in a previous article that the quantum cos- 



mological billiards as a Dirichlet problem implies on the 
one hand the pointwise vanishing of what is commonly 
called the "wavefunction of the universe" as the singular- 
ity is approached, and furthermore that the wavefunction 
is generically complex and oscillating, particularly that it 
can not be analytically extended beyond the singularity 
[16) . In contrast to the standard Wheeler-DeWitt equa- 
tion with a potential, in the quantum cosmological bil- 
liards setting it is consistent to restrict to positive norm 
states, but this requires the wavefunction to be com- 
plex. While complexity of the wavefunction has been 
proposed to be related to the emergence of a directed 
time in canonical gravity [T7HIS], it is also essential for 
a semiclassical wavepacket to have a definite direction of 
propagation. 

The hyperbolic space billiard domains for 13 < 10 pure 
gravity and _D = 11 supergravity, where D denotes the 
number of spacetime dimensions, are non-compact but 
of finite volume in spaces of constant negative curvature. 
The geodesic flow inside these domains is not only ergodic 
and mixing with respect to the Liouville measure, but it 
is furthermore uniformly hyperbolic, having the Anosov 
property [501 HI]. This is the strongest form of chaos, all 
the periodic trajectories are isolated and unstable under 
small variations of the initial conditions. The quanti- 
zation of non-relativistic classically chaotic systems has 
been widely investigated, for such pure reasons as try- 
ing to understand the transition between quantum and 
classical mechanics for chaotic motion, and in particu- 
lar to understand which imprints the nonlinear classical 
systems leave in their respective linear quantum counter- 
parts, see e.g. [32 |23]. This field of research has been 
termed quantum chaos, and billiard systems have served 
as its most popular models. 

Standard semiclassical (or geometric optics) techniques 
are applicable only for classically completely integrable 
systems, and are problematic for ergodic systems due to 
their inconsistency with the isotropic distribution of mo- 
menta in the classical limit, as had already been noted 
in [2T. While the Einstein-Brillouin-Keller (EBK) torus 
quantization of integrable systems is well understood, 
the quantization of classically chaotic systems to date 
requires numerical techniques. Nevertheless, with the 
help of the Gutzwiller-Selberg trace formula (which re- 
duces for integrable systems to EBK quantization) , it has 
been shown that there exist strong correspondences of 
the classical system to its quantized version, relating the 
lengths of the classical periodic orbits to the eigenvalue 
spectrum [211 [55VE7] . Generally, numerically computed 
eigenf unctions of classically chaotic systems show regions 
of increased density along unstable periodic orbits, called 
quantum scars [28j. However, for the bound states inside 
(twice) the billiard domain for D = A pure gravity, quan- 
tum unique ergodicity (QUE) has recently been proven 
[S51 - I5T] . QUE implies that the eigenstates equidistribute 
in the semiclassical limit of high energies, excluding the 
possibility of localization of semiclassical eigenstates onto 
unstable periodic orbits, or of fractal structures. 



In this article, we report investigations concerning the 
behavior of relativistic wavepackets evolving as quantum 
cosmological billiards. Along with analytical considera- 
tions (cf. also |32j), we solve the corresponding Wheeler- 
DeWitt equation numerically, thereby obtaining detailed 
insights into the evolution of semiclassical wavepackets 
and in particular in their long-term behavior. As a pro- 
totype, we study the case oi D — A pure gravity, since the 
higher-dimensional cases are numerically too involved. 
We refer however to [33H37] for results on QUE for higher- 
rank locally symmetric spaces. The main results reported 
in this article are as follows: 

• The relativistic wavepackets follow their corre- 
sponding classical trajectories and then start to de- 
viate due to transversal spreading and reflections 
off the moving billiard walls. 

• The wavefunction tends to zero as the singularity 
is approached. 

• In the flat space description of the cosmologi- 
cal billiards, the energy expectation value of the 
wavepacket decreases due to successive redshifts 
upon being reflected off the moving billiard walls. 

• In the hyperbolic space description, the energy ex- 
pectation value of the wavepacket remains con- 
stant. In agreement with W, we find that a highly 
excited state remains highly excited on its way 
into the singularity. However, the vanishing of the 
spreading wavepacket still suggests a quantum res- 
olution of the singularity. 

• As an aside, we obtain a set of exact solutions for 
the relativistic quantum billiards corresponding to 
the one-dimensional infinite square well with mov- 
ing walls, for the massless as well as for the massive 
case. 

In sectionlTTl we summarize the notation and derivation 
of the Wheeler-DeWitt operator for the setting of quan- 
tum cosmological billiards. In section [Hlj we investigate 
the reflection of a relativistic wavepacket off a moving 
wall using the one-dimensional example of a moving in- 
finite square well potential and provide a set of analytic 
solutions. In section |IV[ we report our results concern- 
ing the quantum cosmological billiards of D = 4 pure 
gravity. These have been obtained through two different 
approaches, namely on the one hand as a quantum bil- 
liards in flat space with moving walls and on the other 
hand in hyperbolic space. We then compare our find- 
ings with the predictions of [TB]. In the appendix, we 
supplement these results with and contrast them to the 
still non-separable but integrable case of the triangular 
billiards in fiat space with all walls kept fixed. 



II. THE WHEELER-DEWITT OPERATOR FOR 
QUANTUM COSMOLOGICAL BILLIARDS 

In the BKL limit, spatial inhomogeneities and matter 
degrees of freedom contained in the original Lagrangian 



description of the gravitational theory are represented 
merely through the infinitely steep billiard walls [71 [16] . 
For the dynamics, one is left with the diagonal part of the 
metric, which can be written for a D = d+1-dimensional 
spacetime in pseudo-Gaussian gauge as 



ds2 ^ -N{TfdT^ + J2 



p-2rM 



(dx^) 



(1) 



/.= ! 



where ^ G {l,...,d}, and where r is any global time 
function parametrizing the foliation into space-like hy- 
persurfaces. We use diagonal coordinates /3'' in the in- 
duced space of the logarithmic scale factors, which we call 
/3-space. In these diagonal coordinates, the Lorentzian 
DeWitt metric reduces to the Lorentzian niinisuperspace 
[14j metric and we can write the corresponding kinetic 
term simply as 



-Ckin - ^'yp-/^"/^' 



(2) 



in terms of the rescaled lapse function n — N/ ^ (see 
e.g. [H]) and with /3'^ = -J^. The spatial volume col- 
lapses to zero at each spatial point at the spacelike singu- 
larity in the past {^/g — >■ 0), where the proper time goes 
towards 0+ . The scale factor Q, of the universe is related 
to the spatial volume ^ through fi = — ^ In ^. We 
choose the lapse for our analysis in flat /3-space such that 
D, = li^ serves as the time coordinate, i.e. the singularity 
is at rj = -|-oo. The remaining /3* with i € {2, . . . , d} serve 
as the parameters of anisotropy. Because of the spatial 
ultra- locality of the BKL limit, the gravitational model 
is reduced to the classical mechanics system of a massless 
relativistic billiard ball which is described by the /3 vari- 
ables and which moves along geodesies with respect to the 
flat DeWitt metric, interspersed by specular reflections 
off the walls. The geodesic parts of the dynamics between 
each two successive reflections are Kasner regimes. For 
each spatial point, there is one such system, and all these 
systems are decoupled from each other. 

With the conjugate momenta tt^ = \'q^vfi'^ ^ where we 
have set the rescaled lapse function ton = 1, we can write 
the free Hamiltonian as ^ = tt^tt^. The billiard walls are 
not all static with respect to Jl-time in the flat Lorentzian 
/3-space, instead some of them move with constant veloc- 
ity. However, by slicing the forward lightcone of /3-space 
in terms of hyperboloids and transforming to hyperbolic 
coordinates (cf. FIG.fl]), we obtain a billiard system with 
motionless walls on a negatively curved space, where we 
adopt the hyperbolic radial coordinate p as the new time 
coordinate. Explicitly, through 



/3' 



PT It^T 



-1 



-M"" 



(3) 



the Wheeler-DeWitt operator can be written in hyper- 
bolic space as 



H = p 



l-d 



d_ 

dp 



,d-i 



d_ 
dp 



iLB • 



(4) 




FIG. 1: Sketch of the cosmological billiards domain as a wedge 
inside the forward lightcone. It is depicted schematically how 
the biUiard wedge is intersected by a spaceHke hyperboloid. 



In the hyperbolic coordinates, the singularity is at 
p = -f(X). We note that the factor ordering is deter- 
mined through the transformation from flat space, i.e. the 
whole Wheeler-DeWitt operator % is simply the Laplace- 
Beltrami operator in hyperbolic coordinates, while 



A 



LB 



v''-'d,{v^-''d,)+v^dl 



(5) 



refers to its restriction to a fixed hyperboloid. We 
choose the generalized upper half-plane with coordinates 
u € 1R''~2, V S IR>o and the Poincare metric dcr^ = 
^^^(dw^ -I- d-u^) as a realization of the hyperbolic space. 
The Wheeler-DeWitt equation reads ?i^ = for the 
wavefunction ^, in flat and respectively in hyperbolic 
space. The infinitely steep potential walls of the bil- 
liard domain are implemented as Dirichlet boundary con- 
ditions on the wavefunction. In hyperbolic space, the 
Wheeler-DeWitt equation separates as 

^{p,u,v)^R{p)F{u,v) (6) 

and the radial part admits the solution [TB] 

R±{p) = p^c±Vfe^-(^)'^°8P (7) 

if the eigenfunctions F solve the eigenvalue equation 

Alb-F(u,«) = -A:2f(u,w) . (8) 

For d > 2, it is apparent from ^ that 'J — > as p — > oo. 
The Dirichlet boundary conditions at the infinite poten- 
tial walls that define the billiard table furthermore imply 
that k^ > ^=^, which implies that the wave function is 
generically complex and oscillating [16) . 



III. RELATIVISTIC WAVEPACKETS IN A 

ONE-DIMENSIONAL INFINITE SQUARE WELL 

WITH A MOVING WALL 

The non-relativistic system of a one-dimensional in- 
finite square well including a massive particle evolving 
according to the Schrodinger equation is one of the most 
elementary quantum mechanical systems. If the potential 
walls are however not static but moving, the situation is 
already much more complicated, since the system is then 
not separable anymore and it is difficult to find exact so- 
lutions. Results for this non-relativistic case have been 
obtained first by Doescher and Rice [55] . 

The relativistic moving-wall case is however a much 
less common object of study, and there appear interest- 
ing subtleties. Since we aim at investigating the quan- 
tum cosmological billiards scenario of a massless quan- 
tum particle evolving according to a Klein- Gordon-like 
Wheeler-DeWitt equation, we at first investigate the one- 
dimensional quantum billiards in an infinite square well 
with one static boundary and one boundary which is 
moving outward at constant velocity v. In order to match 
the notation of the other sections, we use J7 as the time 
and /3 as the space coordinate. Where appropriate, we 
write for subsequent referring the equations in this sec- 
tion such that they are valid for an arbitrary number D 
of spacetime dimensions within A < D < 11. We let 
d = D — 1 denote the number of spatial dimensions. 

The one-dimensional case considered in this section 
then corresponds to d = 2. This system is described 
as the initial/boundary value problem (IBVP) 

*|a^ = ondJ" (9) 

_vI/(r!o,/3) = /(/?) {dt^){no,/3)^9{(3) 

with the infinite square well specified by the domain T = 
[0,L{il)] in terms of the length function L{il) = Lq + 
1/(^1 — flo), where i/ denotes the speed of the receding 
wall, < i^ < 1, and il,/3 e R. The transformation 



/3^/?' 



/3 



L{n) 



(10) 



which is usually applied in such a scenario (see e.g. |39| 
for the Schrodinger case) is of not much help in obtain- 
ing analytic solutions for the relativistic moving-wall sys- 
tem because unlike the first-order Schrodinger case, the 
Klein-Gordon (KG) equation will not separate after the 
transformation (10), although the walls will indeed be 
motionless. 

However, by employing the transformation pi) to hy- 
perbolic coordinates inside the forward lightcone in 1 -I- 1- 
dimensional Minkowski space, we can obtain an infinite 
square well system with static boundaries. Since the 7 
coordinates are constrained to lie on the hyperboloid, one 
of the two coordinates is redundant. We parametrize by 
using coordinates on the one-dimensional analogue of the 



hyperbolic upper half-plane through 



1 



7 



2v 



7 



2v 



(11) 



In the hyperbolic coordinates, we are able to find a full 
set of exact solutions if we require as an intermediate step 
that Lq = vto. We may afterwards freely shift the tip of 
the lightcone in order to obtain solutions for general pa- 
rameters {Lq, to). After the transformation to hyperbolic 
space, the time dependence of the right boundary drops 
out and it depends merely on the constant velocity. 



A(^) = 



1 + U + V1 -ly'^ 



1 -v + ^l-v'^ ■ 
The transformed massless KG equation reads 

pdp {pdp^{p, v)) = vd^ {vdy'^ip, v)) 
and the general solution is given by 



*(/9,-y) = -tpiipv) + ip2{-) 
P 



(12) 



(13) 



(14) 



e.g. 



*n(P,w) 



sin (fc„ ln(u)) exp(-ifc„ ln(p)) (15) 



with 



l^n — 



ln(A(^)) 



(16) 



We note in passing that the massive case corresponding 
to 9q<I> = 5?$ — m?^ can be treated similarly, leading 
to the solutions 

$„(p,w) = Csin(fc„ln(w))[Jifc(mp) + l^fe(mp)] , (17) 

where Jik and Ylfc are Bessel functions of imaginary order 
ifc. 



the Klein-Gordon-li 

(■01, -02) ^ 



The solutions ( 15 ) are normalized to 1 with respect to 
ie scalar product 



ip 



„d-l 



<:\vomdpip2 , 



(18) 



where ip dp (j) 



^dp 



gtp and where dvol = 



denotes the volume element on the (d— 1)- 
dimensional hyperboloid. We can express the solutions 
(15) in flat coordinates as 



^„(f^,/3) 



1 



fc„ln 



X exp 






-ifc„-ln(r!'2 



P' 



(19) 



withr2'(f]) = f7-fio + — ,0<J/<1. The solutions (19) 



are normalized to 1 with respect to the standard form of 
the Klein-Gordon-invariant scalar product. 



('0l|V'2) = 



id-1 



mdn^2 . 



(20) 



Furthermore, the respective norms are preserved in hy- 



perbohc as weh as in flat space. The set of solutions ( 15 ) 



and ( 19 ) seems to be original to us, see however the ap- 



pendix of 40^ for first attempts in the treatment of the 
relativistic one-dimensional moving-wall system. 

For the set of solutions ( 19 ) we can directly show that 



(21) 
A one-dimensional relativistic Gaussian wavepacket 

^{VL,p) = A dpe- " ^ +iW-u,n) (22) 



composed out of positive frequencies uj = \p\ will have a 
norm of one with respect to ( 20 ) if 



A^--[e = P5 + 0fpQC erf(poc) 



(23) 



where erf denotes the error function. The wavepacket 



(22 1 has constant norm (of approximately 1 if c is chosen 



small enough) with respect to ( 20 1 during its evolution in 
the moving wall system, however its absolute width with 
respect to the /3-coordinate will grow with each reflection 
off the moving wall, as the wavepacket will go through 
successive redshifts. If a one-dimensional wave reflects 
off a wall which is moving away at a constant speed, it 
experiences a redshift given by 



1 



./icf 
Zinc 



7^(1 + /^)^ 



1 



1-1^ ' 



(24) 



the redshift is 1 + z 



with"/ — {1 — 1/ ) 2 J e.g. for J/ = |, 
3. In terms of the reduced Hamiltonian H — y^TTg^, 
where 2H = — ttq -|- H^, the corresponding reflection of a 
classical massless particle in one dimension off a moving 
wall implies the relation 



1 + z 



Hr, 



Hi, 



(25) 



which analogously simply states the relation of the energy 
iJjnc of a photon before a bounce from a moving mirror 
to its energy TJ^of afterwards. As a measure for the red- 
shift, we employ the expectation value (H) of the reduced 
Hamiltonian operator, which we call the energy expecta- 
tion value in the following. In order to avoid the square 
root in the reduced Hamiltonian, especially in view of the 
higher-dimensional cases, we adopt the two-component 
notation 41j for ^ by defining 

■0 = </> + X i^o'^ = </> - X (26) 

and adapt it to the massless case to obtain 



H 



(73 + 1(72 



(73 — 1(72 



(27) 



2 2 

where (Ti are the Pauli matrices and where A is the 
Laplace operator in flat space. Upon defining 



^-ii) 



(28) 



equation ( 20 ) is expressed through 

(*i|^2) =2 /"d'^-i/3*l(73*2 
and the energy expectation value through 
(H) =2 I A'^-^p-^^azH-^ . 



(29) 



(30) 



The two-component form (27) implies the expected ex- 
pression 



(H' 



A'^-'prdJM) 



(31) 



for the expectation value of the squared reduced Hamil- 
tonian. Upon reflection off a moving wall in flat space, 
the wavepacket loses energy into the wall, while its KG 
norm is preserved. See FIG. [2]for an illustration of such a 
bounce, calculated by numerically solving the KG equa- 
tion on the one hand in flat space with moving boundary 
conditions, and on the other hand an illustration of the 
corresponding evolution in hyperbolic space with static 
boundary conditions in FIG. [3] In order to observe the 
level excitation distribution for the bound problem in hy- 
perbolic space, we set 



(Hn^) 



-IP 



„rf-l 



dvolV^*9p(ALB0) 



(32) 



and thus the expectation value of the reduced Hamilto- 
nian on the hyperbolic plane in two-component notation 
as 



{Hp) = 2 / dvol^VgiJp* 



(33) 



with 



fr _ rf-2 ^^3 + 1^2 . I 2-rf O"3-lg2 ,o.^ 

iip--p ^ Alb+P ^ • (34) 



For the derivation we have set 



'4j^(t) + x ip"* ^9pV = </- - X 



(35) 



Note furthermore that in the one-dimensional case con- 
sidered in this section, the supremum norm IV'loo of the 
free wavepacket remains constant, in accordance with 
the solution (ItI to the radial part of the d-dimensional 
Wheeler-DeWitt equation for d = 2. 

Using ([29]), we compute position expectation values 
according to 

(/3,) = 2 /"d'^-V*V3A* = i ( A^-^PP^rdn^ (36) 



and analogously in hyperbolic space using ( 18 ) 
mark however that 
nientum operator 



We re- 



unlike the Hamiltonian and the mo- 
the position operator in relativistic 
quantum mechanics generally mixes positive and nega- 
tive frequency components of the wavepacket [41J. We 
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FIG. 2: Plots of |t/'(/3)|^ of a one-dimensional Gaussian wavepacket redshifted upon reflection off a moving wall in flat space. 
The vertical bar in (a)-(c) represents the moving wall. The horizontal scaling in (d) differs from (a)-(c). 
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FIG. 3: Plots of \'ip{v)\^ of a one-dimensional Gaussian wavepacket in an infinite square well in hyperbolic space. 



nevertheless adopt the prescription (36) as the intuitive 



measure for the "center of mass" of the wavepacket and 
employ the usual Klein-Gordon covariant wavepackets 
(cf. e.g. (40) below). For details on localization of rela- 



tivistic particles, we refer the reader to the investigations 
in [321 113 , and to [44^ in the setting of quantum cosmol- 
ogy. 



IV. QUANTUM COSMOLOGICAL BILLIARDS 

The classically strongly chaotic quantum BKL billiards 
corresponding to D = A pure gravity serves as the proto- 
type for studying the relativistic wave chaos inherent also 
in the relevant higher-dimensional theories (see e.g. [?])■ 
If D > 4, the situation becomes quickly unfeasible for 
numerical studies. In the case of D = 4 pure gravity, the 
billiard table is given by a wedge in three-dimensional 
fiat Lorentzian /3-space, cf. FIG. [l] or equivalently by a 
triangle with one of its three edges moving outwards at 
constant speed v as the singularity is approached, see 
FIG. 



4(a) 



It is not known how to solve such a trian- 
gular two-dimensional moving wall system analytically. 
Although by a suitable transformation the moving tri- 
angle walls in two-dimensional flat space can be trans- 
formed to static ones in hyperbolic space, the resulting 
static domain will be of too complicated shape. Thus we 
resort to numerical methods for our investigations. Since 
for the numerical treatment we do not have to rely on 
separability of the Wheeler-DeWitt equation, we studied 



the quantum cosmological billiards in flat as well as in 
hyperbolic space and confirmed that the results match. 



A. Flat space description 

For the case of d = 3, we choose the convenient stan- 
dard notation fi = /?^, /3+ = /3^, /?_ = /3^ for the coor- 
dinates. The Wheeler-DeWitt operator corresponding to 
(B then takes the form 



H = v'^-'d^d, = 









dp-' 



(37) 



and the IBVP of the quantum cosmological billiards is 
defined by 



rH*(f},/3±) = o 

I *|a^ = 
]vI/(r!o,/3±) = /(/3±) 
[(ao*)(f7o,/?±)-.9(/3±) 



in J^ 
on OF 



(38) 



with J- = {/3e R2|^_ >0,^_ <\/3/3+,/3_ < -^{-13+ + 
51)} representing the triangular domain with one face of 
the triangle moving outwards at a constant spee d v = 



dl3\2 



^i=i\ an 



i as displayed in FIG. 



4(a) 



derivation of the wall forms is detailed in [7| |8j. 



The 

Note 

that this triangular domain is actually the sixth part of 
the Mixmaster domain of a Bianchi IX universe 1451 . 
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(a)The cosmological billiard domain for D = 4 pure gravity in 

flat space, here Oi = 32 and Qq = 24. The displayed trajectory 

corresponds to a classical particle with initial conditions 

according to (pTl. The particle starts at the small depicted 

circle and is at first reflected perpendicularly from the receding 

wall at a and then passing again through the starting point. 
Specular reflection does generally only hold in the rest frame of 

the respective wall. After the last depicted bounce at c, the 

classical particle is moving exactly parallel to the /J+ axis and 

will never hit a wall again. 
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(b)The wavepacket's position expectation values according to (J36| 

are centered on the classical trajectory (cf. FIG. [4(a)| | in the 

beginning and then deviate as the wavepacket spreads and 

reflects. This is an enlarged excerpt of FIG. [6(c)] 



FIG. 4: Classical and wavepacket trajectories in the flat space 
description of the quantum cosmological billiards according to 
initial parameters (47 1. 



We rotate the cosmological billiard triangle counter- 
clockwise by a = ^ and perform the time-dependent 
rescalings 



/3; 



p. 



L+m 



p'- 



/3- 



_ \/3f 



(39) 



using the length functions i+(ri) = 2 ^ ^^"^ L_{VL) 



\n 



In the rescaled coordinates (/3:^,/3:L), the quan- 
tum billiards is defined in the static domain given by 
the upper- left half of a square with side length 1, with 
the Wheeler-DeWitt operator modified accordingly. The 
upper-right corner is constrained to be null on the for- 
ward lightcone, i.e. it represents the cusp of the triangle 
on the hyperbolic upper half-plane, cf. the following sub- 
section. 



The Wheeler-DeWitt equation is a second order hy- 
perbolic PDE, thus we require two initial conditions. We 
take a plain Gaussian covariant KG wavepacket and its 
time derivative, chosen such as to restrict to only posi- 
tive frequencies. Since the treatment will in any case be 
limited by the numerical accuracy, the wavepacket will 
fit to any desired accuracy into the box if its width c is 
chosen small enough. A free two-dimensional Gaussian 
wavepacket composed out of positive-frequency modes is 
given by 



*(^,/3±) ^A I dV" 



-4(p- 



-Po)^+ip-(/3- 



-/3o) — iLj(0 — Oq 



)(40) 



with ijj{p) — •\/p+^ +p_^ expressing the relativistic dis- 
persion relation. Using (39), we can write the two initial 
conditions at fig as 



*|no = A 



A^pe 



(41) 



(9a*) k 



-lA 



A^p[u:{p)~p+p+-p^(i4x 



(42) 

where here /3 ^ q^ and where the minus sign in ( |42[ ) 
arises from restricting to positive-frequency modes. The 
condition (41 1 can easily be integrated analytically, re- 
sulting in 



*l 



--.2 



(43) 



however ( 42 ) is not separable in the momentum vari- 



ables and for po ^ has to be evaluated numerically, 
e.g. through Simpson's rule. Requiring the wave packet 
to have unit norm with respect to the KG inner product 



( 20 1 fixes the normalization constant to 

A 1 1 _I ££^72 

A = -c^TT 4e 4 Po X 
2 



(c-2 +p2) J„(^p.2) +p.2^^(^p-.2) -^ ^(44) 



where I^ is the modified Bessel function of the first kind 
and of order k. 

The direction of po determines the direction of prop- 
agation of the wavepacket. Unlike the non-relativistic 
Schrodinger case, the magnitude of po does naturally not 
determine the wavepacket's speed of propagation, which 
is fixed to 1, however it does determine its modulation 
frequency, i.e. its energy, and it effects the rate of its 
transversal spreading. Since we are considering a rela- 
tivistic wavepacket, we have to deal with a dispersion 
relation that states a linear relation between the abso- 
lute momentum and the frequency, unlike the quadratic 
relation in the Schrodinger case. This implies that the 
group speed equals the phase speed and that there is no 



dispersion of the wavepacket in the direction of propaga- 
tion, however there will be a spreading transversally to 
the direction of propagation. Setting higher values for pq 
while keeping c fixed means generating a packet which 
will spread more slowly. 

If a plane wave hits perpendicularly a wall which is 
moving away at a speed z^waib it will experience a red- 
shift 1 -I- z as in (24). As in the one-dimensional case, 
if ^'waii = 5, then 1 + z = 3. This also holds for the 
classical description of the bounce, j^ = 3. However, 
the quantum wavepacket is composed out of a superpo- 
sition of plane waves with slightly non-parallel wave vec- 
tors, representing its transversal localization and fanning. 
Thus, in the wavepacket analogy of a classical perpen- 
dicular bounce, almost all plane wave components are 
in fact not hitting the wall orthogonally, implying that 
the energy will actually be reduced by a factor of less 
than three. For example, for an orthogonal bounce of a 
wavepacket with the initial conditions (47), this results 

in , -^ ° sa 2.89. In the limit of large momenta, we then 

(-H)ret 

recover the classical result of three. 

Specular reflection does not hold for the massless parti- 
cle bouncing off a moving wall. For non-orthogonal such 
bounces, we recall the relations for a classical pointlike 
particle [1] 



Hi, 
H, 



rcf 



sin(6>i.cf) 
sin(6'i„c) 
1 . ,„ 



smlU^ct) 



smi 



(45) 

+ ^inc) , (46) 



where 6'inc and 6'rcf denote the angles of incidence and 
reflection, respectively. For example, if ^inc = f , then 



1.953, while for a wavepacket with \po\ = 5, we 
1.97 if it approaches the moving wall at 



(H) 



obtain ,^,, 

an angle of incidence 6'inc ■ 

For the numerical investigation of the quantum cos- 
mological billiards, we need to specify the parameters for 
the initial data (41 ) which are to be evolved by (38 1. The 



size of the initial billiard domain and the modulation fre- 
quency of the wavepacket effect the necessary grid size 
of the numerical computation. Computational efficiency 
therefore suggests to choose the parameters for the ini- 
tial conditions such that the initial billiard domain is just 



large enough for the initial wavepacket (41 1 with width 
c = 1 to fit in up to numerical accuracy. One possible 
choice, which we adopt for our simulation, is 



f^o = 24 , (3 



+0 



/?- 



Po = Rid , (47) 



where R = -^ 



rotates clockwise by ^. The en- 



Vs 1 

suing wavepacket evolution corresponds to the classical 
path depicted in FIG. 4(a) We have set for our in- 



vestigations IpoI = 5 because the momentum is on the 
one hand large enough such that the Gaussian momen- 
tum distribution is only negligibly indefinite, and on 



the other hand a low initial momentum makes numer- 
ical calculations easier because the necessary grid res- 
olution scales with the modulation frequency squared. 
In the end, we can expect the results of our inves- 
tigations to hold qualitatively also for larger and dif- 
ferently directed initial momentum choices as well as 
larger initial billiard domains. We studied the evolu- 
tion of this two-dimensional relativistic wavepacket nu- 
merically by specifying a two-dimensional grid. The re- 
ported results have been obtained using a grid size of 
5536 X 3200. A three-level convergence test including the 
additional grid sizes 4152 x 2400 and 2768 x 1600 ensured 
us about the convergence of the numerical solution (see 
e.g. [461 El]). The Courant-Friedrichs-Lewy (CFL) con- 
dition (see e.g. [^F.^^P) states that the numerical domain 
of dependence must include the analytical one, implying 
AJ7 < a A/3 with the Courant number cr < 1. We have 
set cr = ^, and thus Ail = 0.0025, and increased Aft 
linearly according to the change of size of the billiard do- 
main. We used fourth order classical Runge-Kutta time 
integration in conjunction with the method of lines (see 
e.g. [50j ) and implemented the spatial derivative opera- 
tors as fourth order finite difference approximations (see 

e.g. m)- 

Snapshots of the ev olution are displayed in FIG. [5J As 
plotted in FIG. |6(b)[ the supremum norm tends to zero 
with growing il, matching the predictions of [16j . In- 
deed, with increasing D, the wavepacket has increasingly 
more space to spread about. We contrast this to the 
case of a static triangular quantum billiards in the ap- 
pendix. We furthermore computed norm, variance and 



expectation values according to (20) and (36), where we 



relied on Lagrange interpolation of the spatial grid for 
the computation of the il derivative. The plot of the 
variance in FIG. 6(e)| suggests growth of the variance 
without bounds as the billiard domain keeps increasing, 
i.e. the wavepacket is spreading out more and more, on 
the one hand because of the transversal fanning, on the 
other hand because of the reflections off the moving wall. 

The Klein-Gordon norm ( [20l) s tays constant within 
bounds of 3 • 10^"*, cf. FIG. |6(a)] For the chosen ini- 
tial parameters, the position (i.e. anisotropy) expecta- 
tion values computed according to ( 36 ) increase towards 



the cusp as the singularity is approached, as plotted in 
FIG. |6(d7] However, this can be seen to happen because 
the corner of the triangle which corresponds to the cusp 
on the hyperbolic upper half-plane is on the lightcone and 
thus the part of the wavepacket moving into this direc- 
tion will only slowly start to reflect from the moving wall 
and also spread more slowly, thereby gaining in relative 
weight in the anisotropy expectation value computation. 
We furthermore observe by comparison of FIG. |4(b")| with 
FIG. 4(a) that the wavepacket's trajectory in terms of 



its anisotropy expectation values ( 36 1 follows at first the 



classical path and then deviates due to the spreading of 
the wavepacket towards the walls. 

The energy expectation value and variance are given 



n=24.0006, <v|\(()=0.999941 

<p+)=7.99968, <P_>=4.5, 0^=0.994523 

supM^=0.0315142, <H)=5.15674, (AHf =0.477618 




12=32, <\((|\|f)=0.999958 

(P+ >=1 1 .5525, <P_ >=1 0.6522, o-=2.77641 

sup|\(<p=0.0593349, (H>=3.78938, (AH)^=3.0213 
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(a)Initial circular Gaussian wavepacket 
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(b)The wavepacket near point a (cf. FIG.[4(a)]l. 
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(c)The wavepacket near point c (cf. FIG.[4(a)||. 
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(d)The wavepacket is moving towards the cusp while spreading 
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(e)At large il, the wavepacket has spread over the billiard domain. (f)Magnification of the cusp in FIG. [5(e)] Since the right corner is 

on the lightcone, the wavepacket can never reach it. 
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FIG. 5: Snapshots of wavepacket evolution in the quantum cosmological billiards of D = 4 pure gravity in flat space. 



through (30 1. From the data displayed m FIG. |6(f)| one 
observes the decrease in energy due to the wavepacket 
getting redshifted successively. In the following, we com- 
pare our results to the semiclassical considerations of [1] 
for the Mixmaster universe. Although the billiard do- 
main considered here is only the sixth part of the Mix- 
master billiard domain, we can nevertheless relate the 
respective results through the observation that the Mix- 
master domain can be folded three times along its three 
static heights to match the billiard domain considered 
here. Though a comparison should reveal a difference in 



the rate of decrease of the supremum norm and in the rate 
of increase of the position expectation values and vari- 
ance, we can expect the results about the rate of energy 
loss to be the same in both cases. Through an adiabatic 
approximation, Misner examined in T on the one hand 
the bound states of the billiard triangle for a given instant 
of time J7 and on the other hand the classical trajectories 
of massless particles and the energy loss involved from 
reflections off moving walls. From this analysis it turns 
out that the product of energy expectation value and Q is 
roughly invariant, and thus that a semiclassical, highly 
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(a)The norm of the wavepacket during its evolution according to 

( |20| . The parabolic envelope is due to the enlarging grid spacing, 

the oscillations mainly stem from interpolation artefacts in the 

computation of the Q derivative. 
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(c)Extended parametric plot of the of the position expectation (d)Long-term behavior of the expectation values of the anisotropy 
values according to ||36|, cf. FIG.[4(b)| parameters 
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(f) Energy expectation value of the wavepacket 
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FIG. 6: Data of the wavepacket evolution with initial values (47 1 in the quantum billiards of D = 4 pure gravity in flat space 
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excited state has to remain highly excited all the way 
into the singularity. Comparing the results of [4] with 
our results displayed in FIG. 6(h) , we indeed find that 
in the long term, the product {H)Q stops increasing and 
tends to remain constant, in fact more and more as the 
localization of the wavepacket is lost. Our results thus 
suggest that the results of the semiclassical analysis of [4| 
are smoothed out via an infinity of classical trajectories 
realized at the same time, retaining the constancy of the 
product {H)D,. In terms of wavepackets, this in turn sug- 
gests that the constancy of (H) il amounts to a complete 
loss of localization in the sense that the wavefunction is 
evenly distributed across the billiard domain. 

In the following section, we report the results of our 
investigations of the quantum cosmological billiards in 
hyperbolic space. Since the billiard walls are static there, 
an adiabatic approximation will not be necessary in the 
first place. 



B. Hyperbolic space description 

We induce a coordinate transformation to hyperbolic 
space with p as the radial coordinate by setting 



7" 



is then expressed via 



p-'P"- 



(48) 



7 = 



7 



ird domain of D = 


4 pure gravity 


(u-^r+v^ + i 




VSv 




(u-'^r+v'-i 


(49) 


V^v 


\-u 





on the upper half-plane (UHP). The now static billiard 
walls are accordingly specified by 



u = 



1 



.^^1, 



(50) 



as opposed to the moving domain walls in /3-space (cf. 
FIGS. [4(a)] and [7(a)| ). The biUiard domain on the UHP 
as specified by (|50[) coincides with the fundamental do- 
main of the group PGL2(Z), which is half the fundamen- 
tal domain of the modular group PSL2(Z). For numeri- 
cal investigations of the purely discrete spectrum of odd 
Maass waveforms for PSL2(Z) and PSL2(Z[i]), we refer to 
[521 - I55] . The initial parameters pTJ) of the previous sec- 
tion translate via the relations ( |49| ) to hyperbolic space 
and are denoted by {pq,uq,vq) in the following. 

The initial wavepacket, which is a modulated Gaussian 
in flat space, is expressed on the hyperbolic plane through 



^{po,u,v)^A dV 



"Irip-Po? 



^Pof{p+,P-,u,v) 



(51) 



with f{p+,p_,u,v) given by 

f{p+,p-,u,v) =p+{'^^{u,v) -7^(uo,wo)) 



(52) 



The evolution equation for the quantum cosmological 
billiards is determined by the Wheeler-DeWitt equation 
with the operator Q specified to d = 3, namely 



'o7 



2p 



d_ 
dp 



A 



LB 



'^{p,u,v) = 



Here, Alb denotes the Laplace-Beltrami operator 



^LB 



du 






(53) 



(54) 



on the UHP. As in the previous section, we translate the 
potential walls into Dirichlet Boundary conditions and 
study the evolution of the wavepacket numerically, where 
we implement the differential operators through a finite- 
differencing scheme using the method of lines, of fourth 
order in the space-like variables (u, v) as well as in the 
time-like variable p through fourth order Runge-Kutta 
time integration. 

In order to apply the numerical method to the hy- 
perbolic description of the quantum cosmological bil- 
liards, the second order Wheeler-DeWitt equation (53) 



is split into a system of two first-order PDEs by defining 
^p{p,x,y) = p'^£-^{p,x,y), leading to the IBVP 



P 
*(p,u,w) : 

^*(po,u,f) 



= Alb*(/5, u,w) in J" 
in T 

'ip{p,u,v) — on dF 
:/(u,w), ■^p{pQ,u,v)= g{u,v) 



2j^*(p,u,t-) = *p(p,u,w) 



(55) 



In addition to ( |51| , a second initial condition is re- 
quired. The second one, pQ^{dp'^){pQ,u,v), is given by 

/CO ^ 

A'pc-'^^P-P^<^"f{p+,p.,u,v)x 
-OO 

X [p+-i''{u,v) + p--i'^{u,v) - uj{p+,p-)-i^{u,v)\ . (56) 

For the efficient implementation of the numerical 
method, we transform the hyperbolic billiard domain 
(cf. FIG. 7(a) ) into a rectangle and then compactify 
it. Our results have been obtained using a grid size 
of 1600 X 3920. They have been subjected to a three- 
level convergence test including the additional grid sizes 
1200 X 2940 and 800 x 1960 in order to check the vafidity 
of the results. 

These results show that, similar to the fiat space de- 
scription, the wavefunction tends to zero pointwise with 
increasing p in the hyperbolic description of the quan- 
tum cosmological billiards, cf. FIG. |9(b)[ This sug- 
gests the vanishing of the wavefunction at the singularity, 
matching our predictions of |16j, see also equation ([7|. 
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FIG. 7: Classical and wavepacket trajectories in the hyper- 
bolic space description of the quantum cosmological billiards 
according to initial parameters (471. 
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(b) Parametric plot of the 

trajectory in terms of t he 

expectation values, cf. FIG. [7(a)] 



The position expectation values according to ( |36[ ) of the 
wavepacket 's evolution in the billiard domain on the hy- 
perbolic UHP are displayed in FIGS. [9(c)1 and [9(d)| The 
parametric plot of the position expectation values dis- 
played in FIG. [7(b")1 shows that the wavepacket follows at 
first the classical trajectory shown in FIG. 7(a) and then 
starts to deviate from it as it is spreading out. After the 
last reflection at point c, the classical path depicted in 
FIG. 7(a)| is following a straight "vertical" line towards 
the cusp. Since on the UHP every geodesic with respect 
to the Poincare metric is either a half-circle or a verti- 
cal straight line, all the neighboring rays representing the 
wavepacket are given by half-circles. However, all these 
rays are subject to reflections from the billiard walls and 
will only reach a finite height with respect to the v co- 
ordinate. This is reflected in the wavepacket's trajectory 
displayed in FIG. |7(b)[ which has a turning point at a 
certain value of v, contrary to the classical path. Con- 
tour plots of the evolution of the wavepacket are shown 



in FIG. [8] at different instants of p, chosen such as to 
display similar situations as in FIG. [5J 

The energy expectation value and its variance de- 
scribed in the hyperbolic description by ( 33 1 converge 



with increasing numerical resolution to a constant value 
throughout the wavepacket's evolution, cf. FIGS. [9(f)[ 



and 9(g) We pick up where we have left at the end of 
the preceding section and compare our results with the 
semiclassical considerations of [J. In the hyperbolic de- 
scription of the cosmological billiards, the billiard walls 
are static and thus we do not have to rely on an adiabatic 
approximation, contrary to the description in flat space of 
the previous section. Since the energy expectation value 
and its variance remain constant, we directly see that a 
highly excited semiclassical state stays highly excited on 
its way "into" the singularity as p — > oo. On the other 
hand, we observe that the localization of the semiclas- 
sical wavepacket is lost due to the spreading across the 
billiard domain, and particularly that the wavefunction 
tends to zero as the singularity is approached. In [1] it 
was deduced that the quantum state remains classical all 
the way into the singularity. In the current framework 
however, we propose to rather take the vanishing of the 
wavefunction due to its spreading (and thereby its loss of 
localization) than the excitation in terms of energy levels 
as an indicator for a quantum resolution of the singular- 
ity. 

We conclude this section by remarking that the 
wavepacket evolution in hyperbolic space can be mapped 
back to flat space in order to double-check the correct- 
ness of the computation. The set of hyperboloids which 
make up the discretization of the forward lightcone has to 
be sliced into equal-i7-time hypersurfaces, such that for 
every such "horizontal" slice, there is a one-dimensional 
contribution from each intersecting hyperboloid. Using 
interpolation, one can then confirm that the evolution un- 
der (53) with the initial conditions (51) indeed matches 



the evolution under ( 38 ) with the initial conditions ( 47 ) 



that has been reported in the previous section. 



V. SUMMARY AND DISCUSSION 

Pure gravity and supergravity in four to eleven space- 
time dimensions admit a description in terms of cosmo- 
logical billiards in the BKL limit towards a spacelike sin- 
gularity. These singularities of the classical theory are 
expected to be resolved through quantum effects. The 
classical billiard motion of the BKL limit is classified 
as strongly chaotic for theories such as four-dimensional 
pure gravity and eleven-dimensional supergravity, and 
exact solutions for the quantization of these classically 
chaotic billiard systems are not available. With this ar- 
ticle, we supplement and extend our previous results [TB] 
by reporting investigations regarding the evolution of rel- 
ativistic wavepackets towards the singularity as the quan- 
tum cosmological billiards associated with D = A pure 
gravity. In agreement with our previous results |16j . we 
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FIG. 8: Snapshots of wavepacket evolution in the quantum cosmological billiards for _D = 4 pure gravity in hyperbolic space. 



found that the wavefunction of the universe decreases 
pointwise towards zero as the singularity is approached, 
while its norm is preserved. 

The initially localized wavepacket evolves according to 
the Wheeler-DeWitt equation, which reduces to a mass- 
less Klein-Gordon equation for the finite-dimensional 
cosmological billiards. The subtleties of relativistic 
wavepackets in domains with moving walls have been 
elucidated first in the one-dimensional case in section 
|III[ and an exact set of solutions has been derived, for 
the massless as well as for the massive relativistic Klein- 



Gordon quantum particle inside the infinite square well 
with moving walls. Since analytic solutions for the quan- 
tum cosmological billiards associated with D = A pure 
gravity are however not known and out of reach, we in- 
vestigated the relativistic wavepacket evolution using nu- 
merical methods in section HVl Two different realizations 
have been employed. The first one is in flat space with 
moving billiard walls, while the second one is in hyper- 
bolic space with static walls. The solutions in the two 
realizations for corresponding choices of initial parame- 
ters can be mapped into each other, and in both cases 
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FIG. 9: Data for the wavepacket evolution according to initial values (471 in the quantum cosmological billiards for D = 4 pure 
gravity in hyperbolic space. 



tend towards zero with preserved norm as the singularity 
is approached. 

In the flat space description of the cosmological bil- 
liards, the position (i.e. anisotropy) expectation values 



follow at first the corresponding classical trajectory and 
then start to deviate as the wavepacket spreads across the 
billiard domain with increasing variance due to transver- 
sal fanning and reflections off the moving walls. Upon 
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such reflection off a moving wall, the wavepacket loses 
energy into the wall and gets redshifted. As in the flat 
space description, the position expectation values of the 
wavepacket in the hyperbolic space description of the 
quantum billiards follow at first the classical path be- 
fore deviating due to the spreading. A classical particle 
can travel on a geodesic all the way into the cusp of the 
fundamental domain at infinity, whereas the wavepacket 
trajectory in terms of its position expectation values nec- 
essarily has a turning point. The corresponding sections 
of this article have each been supplemented with a set 
of snapshots of the evolution of the wavepacket. These 
results can be contrasted to the quantum billiards of rel- 
ativistic wavepackets in a static triangle, reported in the 
appendix. 

We furthermore observe that the product of energy ex- 
pectation value and time tends to a constant value in the 
flat space description, in agreement with the results of the 
semiclassical analysis of ^ which imply that a highly ex- 
cited semiclassical state remains highly excited and thus 
classical all the way into the singularity. However, the 
preserved excitation becomes more obvious in the hy- 
perbolic space description of the cosmological billiards, 
where the billiard domain of 13 = 4 pure gravity happens 
to coincide with the fundamental domain of the group 
PGL2(Z). There, the billiard walls are static, making an 
adiabatic approximation redundant, and we found that 
the energy expectation value and its variance are pre- 
served. Nevertheless, regardless of how highly excited in 
terms of eigenstates, the loss of localization of the initial 
wavepacket as it spreads across the billiard domain, and 
particularly the vanishing of the wavefunction suggest 
the interpretation that the singularity is resolved via an 
effective disappearance of spacetime near the singularity, 
making it effectively unreachable, possibly in line with a 
"deemergence" in terms of a Kac-Moody algebra coset 
model description [56 . In this light, we thus propose 
that the vanishing of the wavefunction and its spreading 
across the whole billiard domain may be considered as 
more appropriate indicators for a quantum resolution of 
the classical singularity than the excitation of the quan- 
tum state in terms of energy levels. 

It should be remarked at this point that the quantum 
mechanical avoidance of classical cosmological singulari- 
ties has also been discussed in terms of approaches that 
are different from the one employed here. In some of 
the models recently investigated in [S7j , vanishing of the 
wave function near the classical singularity is found, too. 
In general, the resolution of the singularity that is sug- 
gested by the vanishing of the wavefunction is to be con- 
trasted to other mechanisms such as the Hartle-Hawking 
no-boundary proposal 58J or the cosmic bounces of loop 
quantum cosmology (see [59^ for a recent review), since 
the latter two require continuation of the wave packet 
into the singularity and beyond. It should also be noted 
that there exist approaches that circumvent the assump- 
tions of the singularity theorems of Hawking and Pen- 
rose in the first place, for example through a violation of 



the null energy condition from a ghost condensate |60j . 
Such a mechanism can produce a smooth, non-singular 
bounce where an ekpyrotic contracting phase is joined 
with the standard expanding phase [5THB^ , thereby pro- 
viding new solutions for cosmological issues such as the 
flatness problem. 

The setting considered in this article only provides a 
very first step. Albeit very useful, the cosmological bil- 
liards scenario with its infinite potential walls represents 
only crudely the spatial inhomogeneities and possible 
matter degrees of freedom of the theory. Since the shape 
of the billiard table points to an underlying structure in 
terms of Lorentzian Kac-Moody algebras, the quantum 
state of the universe would have to be generalized ac- 
cording to the "small tension" expansion of [BS] as the 
ensuing first step away from the BKL limit of the classi- 
cal theory. We conclude this article by emphasizing that 
it is conjectured that in the deep quantum regime, space- 
time based quantum field theory is replaced by a coset 
model built from the respective Kac-Moody algebra that 
defines the quantum cosmological billiards [Ml IM] • 
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Appendix A: The triangular quantum billiards with 
static domain walls 

In this appendix, we study the triangular quantum bil- 
liards with the same parameters as in the quantum cos- 
mological billiards before, but this time with every single 
one of the billiard walls remaining static from the out- 
set. In this case, the corresponding classical billiards is 
completely integrable. Note that for triangular billiards 
with static walls, only three cases are classically com- 
pletely integrable, namely the equilateral triangle, the 
half-equilateral triangle, or the half-square. These are the 
cases that feature a second constant of the motion, and 
the invariant surface of the phase flow is given by a two- 
dimensional torus in the four-dimensional phase space. 
Despite the non-separability of the triangular billiards, 
analytic solutions are indeed known for the quantum bil- 
liards that correspond to these classically integrable cases 
[66VI69J . For a right triangle with rational angles other 
than the two cases mentioned above, the classical bil- 
liards is only pseudointegrable, and although its trajec- 
tories are still conflned to a two-dimensional surface in 
phase space, this surface is not a torus but has higher 
genus [70j . and analytic solutions for the corresponding 
quantum billiards are not known. We also note that for a 



16 



£i=9884, <>|/|>|»>=0.999961, |>]/r=0.999878 

<|3+>=9.61572, <|3_)=3.54416, 0^=30.2999 

sup|\|/M.00551092, <H>=5. 14822, (AH)2=0.488108 




(a)At f! = 9884, the initially localized wavepacket has spread over 
the static triangular billiard domain. 
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FIG. 10: The quantum billiards in the nonintegrable static triangle that corresponds to the initial cosniological billiard domain 
from the previous section with each billiard wall kept fixed. The norm (vl/l^) ~ |'I'|^ ~ 1 during the whole evolution, oscillating 
only within bounds of 5 ■ 10^®. The grid size was 800 x 800 data points. 



right triangle with non-rational acute angles, there exists the classical motion is ergodic and fills the whole phase 
no other constant of motion besides the Haniiltonian and space. 
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We set the parameters as in (47) and investigate the 



wave packet propagation numerically. As in the calcu- 
lations reported in the previous sections, a three-level 
convergence test was applied on the numerical computa- 
tions. A snapshot of the long-term behavior is plotted 
in FIG. 10 The Klein-Gordon norm according to (20) 
remains constant within bounds of 5 • 10~^, and as ex- 



pected, the variance now has an upper and the supremum 
norm a lower bound. The position expectation values ac- 



cording to ( 36 1 are oscillating around the center of mass 



of the triangle. The energy expectation value according 
to ( 30 1 converges with increasing resolution to a constant 



value, and so does its variance. 
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